Topological states and braiding statistics using quantum circuits 
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Using superconducting quantum circuits, we propose an approach to construct a Kitaev lattice, 
i.e., an anisotropic spin model on a honeycomb lattice with three types of nearest-neighbor interac- 
tions. We study two particular cases to demonstrate topological states (i.e., the vortex and bond 
states) and show how the braiding statistics can be revealed. Our approach provides an experimen- 
tally realizable many-body system for demonstrating exotic properties of topological phases. 



Topological quantum systems are currently attracting 
considerable interest because of their fundamental im- 
portance and potential applications in quantum comput- 
ing 1]. They exhibit topological phases of matter that 
are insensitive to local perturbations. This exotic prop- 
erty makes them appealing for fault-tolerant quantum 
computing. 

Topological quantum computing can be implemented 
by storing and processing quantum information with 
anyons which are neither bosons nor fermions, but 
obey anyonic braiding statistics. To achieve this, it is 
a centrally important issue to experimentally realize a 
topological quantum system. Recently, topological quan- 
tum computing using non-Abelian anyons in a fractional 
quantum Hall system was proposed [2]; however, it has 
not been experimentally verified if the observed fractional 
quantum Hall states are the desired topological states 
for quantum computing. Indeed, up to now, anyons 
have never been directly observed in experiments. Thus, 
in order to guarantee anyonic topological phases in a 
many-body quantum system, one could design artificial, 
but exactly solvable lattice models that have such de- 
sired topological phases. The most important example 
is the Kitaev model on a honeycomb lattice which 
is an anisotropic spin model with three types of nearest- 
neighbor interactions. Most strikingly, depending on the 
bond parameters, this model 's'l supports both Abelian 
and non-Abelian anyons. Therefore, if this topological 
model could be realized in experiments, it would provide 
exciting opportunities for experimentally demonstrating 
anyons and implementing topological quantum comput- 
ing. 

Proposals were made for implementing the Kitaev 
model using atomic optical lattices 0, Also, it was 
shown that optical lattices can be used to demon- 
strate the anyonic braiding statistics. However, its imple- 
mentation using an optical lattice requires extremely low 
temperatures 4] that are beyond current experimental 
capabilities. Moreover, to demonstrate anyonic braiding, 
the optical lattice needs to perform rotations on single 
atoms at selected sites using external laser beams. This 
is difficult because the atomic spacing is of the order of 



laser wavelength and the diffraction limit takes effect. 

Here we propose an approach to construct a Kitaev 
lattice using superconducting quantum circuits 0]- A 
charge qubit is placed at each site and the three types of 
nearest-neighbor interactions are achieved using different 
circuit elements. This circuit can demonstrate topologi- 
cal states (i.e., the vortex and bond states) and show how 
the braiding statistics can be revealed. Here, the quan- 
tum circuits behave like artificial spins [tI and the Kitaev 
model can be constructed using the experimentally avail- 
able technologies for superconducting qubits. Moreover, 
these charge qubits are tunable via external fields, mak- 
ing the required single-qubit rotations implcmentable for 
demonstrating the braiding statistics of the topological 
states. Our proposal provides an experimentally realiz- 
able many-body system for demonstrating exotic proper- 
ties of the topological phases. 

Kitaev model. — We use superconducting quantum cir- 
cuits to construct an artificial Kitaev lattice (see Fig.[T]). 
The building block of the lattice consists of four charge 
qubits [Fig. [ija)] , each placed at a lattice site of a hexag- 
onal lattice [Fig. [Hb)]. Here, each charge qubit is a 
Cooper-pair box connected to a superconducting ring by 
two identical Josephson junctions. Each qubit is con- 
trolled by both the magnetic flux <i>i piercing the SQUID 
loop and the voltage Vi applied via the gate capaci- 
tance Cg. Moreover, as shown in Fig. [T]Ja), the nearest- 
neighboring charge qubits in the three different link di- 
rections are coupled by different circuit elements: 

(i) X-type bond. The two qubits [denoted by 1 and 2 in 
Fig. [TJa)] in the x-link direction are coupled by a mutual 
inductance M. This inductive coupling is given by Q 
Kx{l,2) = MI1I2, where the circulating supercurrent 
li in the SQUID loop of the ith charge qubit is li — 
/csin(7r<i>i/<i>o)cos(/3j, with Ic = 2ttEj/(^q, <i>o = h/2e 
(the flux quantum) and Lpi being the average phase drop 
across the two Josephson junctions. 

(ii) Y-type bond. The two qubits 1 and 3 in 
the w-link direction are coupled via an LC oscilla- 
tor [9| and the inter-qubit coupling is Ky{l,3) = 
-~A(Eji{'^i)Ej3{<^3) sin ipi sin ip3, where Ejii'Pi) = 
2£:jcos(7r$,/$o), and^ = L[7^CE(Cg + C„)/A$o]^ with 
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FIG. 1: (Color online) (a) Schematic diagram of the basic 
building block of a Kitaev lattice, consisting of four super- 
conducting charge qubits (labelled 1 to 4): (i) The qubits 1 
and 2 are inductively coupled via a mutual inductance M; 

(ii) the qubits 1 and 3 are coupled via an LC oscillator; and 

(iii) the qubits 1 and 4 are capacitively coupled via a mu- 
tual capacitance Cm- Inset: These three types of inter-qubit 
couplings are denoted as X, and Y and Z bonds. Here each 
charge qubit consists of a Cooper-pair box (green dot) which 
is linked to a superconducting ring, via two identical Joseph- 
son junctions (each with coupling energy Ej and capacitance 
Cj), to form a SQUID loop. Also, each qubit is controlled 
by both a voltage Vi (applied to the qubit via the gate ca- 
pacitance Cg) and a magnetic flux $i (piercing the SQUID 
loop), (b) A subset of the Kitaev lattice (honeycomb lattice) 
constructed by repeating the building block in (a), where a 
charge qubit is placed at each site. Also, the plaquette op- 
erator Wp — (jf a|a|(jf (Tg(T| is shown for a given plaquette 
(hexagon) p. 



Cs — 2Cj + C„ 
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Note that 



^ oc (Cg + C„i)^ in the presence of the mutual capac- 
itance Cm that connects qubit 1 (3) with its nearest- 
neighboring qubit in the vertical (z-link) direction. Usu- 
ally, Cm 3> Cg , so the coupling between qubits 1 and 3 is 
now much increased, compared to the usual case without 
the mutual inductance 9] , where ^ oc C^ . 

(iii) Z-type bond. The two qubits 1 and 4 in the z-link 
direction are coupled via the mutual capacitance Cm 

K^{1,4) = Em{ni - ngi)(n4 - 71^4), where Em = 



4e^Cm/A, ngi = CgVi/2e, and Ui = —id/dipi is a number 
operator of the Cooper pairs in the ith box. The mutual 
capacitance Cm modifies the charging energy 2e^/CE of 
an isolated charge qubit to become Ec = 2e^Cs/A. 

When the electrostatic and the Josephson couphng en- 
ergies of each charge qubit are included, the Hamiltonian 



of the lattice is H = J2i[-^c{f 



f -Ej,{<^i)cos 



Ex-link ^^(j'fc) + Ey-li„k^y(j'^) + Ez-link^^O''^)- 

For charge qubits, E,. ^ Ej. When the gate voltage 
Vi is near the optimal point e/Cg, i.e., Ugi ~ i, only two 
charge states and corresponding to zero and one 
extra Cooper pairs in the box, are important for each 
qubit. In the spin-i representation, based on the charge 
states |0)i = |t)i and = one has rii — i(l — a^), 
cos (fi — htyf, and sin(/3j — 
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Here we consider the 



simple case with Ugi ~ Ug (i.e., Vi = Vg) and $i 



for all qubits. 
reduced to 



The Hamiltonian of the system is then 
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X— link y — link z — link 

+ J2{h.at + h,a:). (1) 

i 

This is the Kitaev model on a honeycomb lattice, in the 
presence of a "magnetic field" with z- and x-components. 
Here ft,^ = {Ec+^Em)ing-^), and = -^Ej{<^e), with 
Ej{^e) = 2£',/ cos(7r$e/'I'o)- The bond parameters are 

J, = iM/2sin2(^$,/$o) >0,Jy = -^Eji^e)? < 0, 

and Jz — J Em > 0. The magnetic field can be used to 
achieve single-qubit rotations for demonstrating anyonic 
braiding statistics. 

Topological excitations. — Below we focus on two par- 
ticular cases to show the properties of topological excita- 
tions in the Kitaev model ([1]): 

(i) Kitaev lattice in a weak magnetic field. We first con- 
sider the case with |/iz|, l/ial ^ Jx,\Jy\,Jz, and choose 
V = ^i{hz(Tf + hxuf ) as the perturbation. Using per- 
turbation theory in the Green function formalism [3| , one 
can construct an effective Hamiltonian HcS acting on the 
vortex-free sector: Hcs = -(S/i^/Ae^) Ez_iink '^I'^fe ~ 
{2h1/Aex) Ex-link '^j'^k^ whcrc Sz(x) is the excitation en- 
ergy of the state a^'^^^g), and \g) is the eigenstate of all 
plaquette operators Wp [see Fig. Hfb)] corresponding to 
the eigenvalue Wp — 1. Here the effective Hamiltonian 
-ffcff is only contributed by the second-order term because 
both the first- and third-order terms are zero. With the 
zeroth-order term (unperturbed Hamiltonian) included, 
the Hamiltonian of the system can be equivalently 
written as = 4 Ex-iink ^^J^f + Jy J2y -link '^j'^l + 

Ez-iink "^I'-'^fc' where the bond parameters Jz and Jx 
are renormalized to J'^ — Jz ~ 2h1/Aez, and = 
Jx - 2hl/Aex. 

To show the topological excitations, we focus on the 
Abelian case Q with J'^ , | | . The dominant part of 
H, i.e., Hq = J'^ Ez-iink ^j^L favors a highly degenerate 
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FIG. 2: (Color online) Schematic diagram of the procedures 
for braiding topological excitations, (a) The operations Uh 
and Uv for creating different types of topological excitations, 
which are achieved by successively applying spin-pair opera- 
tors at z-bonds along the horizontal [Ph) and vertical (Pi,) 
paths, (b) A combined operation U^^U^^UhUv for both 
braiding the topological excitations created in (a) and fus- 
ing them to the vacuum, (c) The operations Uh and Uv for 
creating the same types of topological excitations, which are 
achieved by successively applying spin-pair operators at z- 
bonds along Ph and Pv. 

ground state with each pair of spins in each z-link aligned 
opposite to each other (|ti) or |it))- While preserving 
the ground-state subspace, the effective Hamiltonian can 
now be written, up to fourth order, as 

H,s ^ ~J,eJ2Wp, (2) 

p 

where Jcs — J'^Jy /WJ'l. Clearly, the ground state of 
the system preserves Wp — 1, i.e., tVpl^) = \g), for all 
plaquettes p. A pair of vortex excitations are created 
when the system changes with Wp — I ^ —1 for two 
neighboring plaquettes. This can be achievedjDy acting a 
spin-pair operator on the ground state \g): \Zi) = a^lg), 



and \Yi) = trf with a^ = afli and af = afaf. Here 
the two operators erf (erf) and li (af) act on the ground 
state \g) at the bottom and top sites of the ith z-link, 
respectively. This pair of vortex excitations are Abelian 
anyons 0, Q , which have an energy gap Ae = 4 Jeff above 
the ground state. 

(ii) Kitaev lattice with dominant z -bonds. We choose 
hz — 0, which corresponds to the case with each charge 
qubit working at the optimal point ng = ^. More- 
over, we choose suitable circuit parameters to have Jz S> 
\hx\, JxAJy\j with I /la; I comparable to Jx,\Jy\- Here we 
also use perturbation theory in the Green function for- 
malism to derive the effective Hamiltonian. Up to the 
second-order term, the effective Hamiltonian can be writ- 
ten as 

HeS — -Keff ^ CTjO"^, (3) 
z — link 

with Kctf — h^/Jz- The spin-pair operator crjcrf at a 
z-bond is commutative with the unperturbed zero-order 
Hamiltonian Hq = JzJ2z-iink'^j'^k' anticommuta- 
tive to the four plaquette operators Wp connected to this 
z-bond. Similar to Wp, the pair operator crjcr^ also has 
two eigenvalues Pz — ±1. Now the ground state \g) of 
the system preserves Pz = 1, i.e., crjcr^ \g) = \g), for all 
z-bonds. When the pair operators af and af are sep- 
arately applied to the ground state at the ith z-bond, 
the excited states \Zi) = af\g), and \Yi) = aflg) are two 
types of bond states; each corresponding to the change 
= 1 ^ — 1 at the zth z-bond and having an energy 
gap 2Kcf{ above the ground state. 

The braiding of topological excitations. — As shown in 
6], an e-type vortex looping around another e vortex 
does not produce a sign change to the wave function, 
while an e vortex looping around an m-type vortex does. 
This indicates anyonic statistics between the e and m 
vortex states. Here we show an alternative procedure 
for braiding a topological excitation with another, which 
applies to both the vortex and bond states. 

First, successively apply spin-pair operations = 
erf erf to the ground state \g) at three z-bonds in the verti- 
cal path Py [see Fig.l^Ta)]. In the case of Hamiltonian 
this creates a pair of e vortex states and move a vortex 
along the vertical path Py , while a bond state is generated 
in the case of Hamiltonian ([3]). Then, we successively 
apply the operations rrf = afli at four z-bonds along 
the horizontal path Ph [see also Fig. EJa)]. After these 
operations, the state of the system is UhUv\g), where 
Uh = (T|(T|(5-|o-f , and Uv = o'sctIo'i correspond to the 
operations along the horizontal and vertical paths. Fur- 
thermore, successively apply U~-^ and Uf^^ to the system, 
so as to fuse [ll,0| the excitations to the vacuum (i.e., the 
ground state) [see Fig.[2fb)]. Now, the final state of the 
system is l^*/) = U^^U~^UhUy\g). Because the paths 
and Ph intersect at a lattice point, where and (tI an- 
ticonimute, one has UhUy = —UyUh. Therefore, the final 
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state becomes \^ f) = —\g)- In contrast, for two topolog- 
ical excitations of the same type, when similar operations 
are applied, the paths Py and Ph do not intersect at a 
lattice point [see Fig. ^c)]. Thus, UhUy = UyUh, and 
I*/) = Id}' yielding no sign change to the ground-state 
wave function. 

The braiding statistics of topological excitations can 
be revealed by means of Ramsey- type interference 0, [l^ . 
To achieve this, we keep the same Uy as above, but use 
Uh = where (CTf)5 = (crf)^/,, 

i.e., each erf is replaced by half of the rotation. In the 
braiding case shown in Fig. [2l^a), ^KctI)^ — i{^2)~^^2 
at the crossing point of paths Ph and Py. Thus, f) — 
U^^U-^UhUy\g) = (a|)-^[z(af)-5]|.g) = zjZa), simi- 
lar to the case with an e vortex looping around a su- 
perposition state of an m vortex and the vacuum 
However, in the case without braiding [see Fig. EJb)], 
1^-^) = C/-^f/-iJ7?,[/„|3) = \g). Therefore, the braiding 
of topological excitations^ can be distinguished by verify- 
ing if an excited state [^2) occurs at the crossing point 
of paths Pfi and Py. 

Discussion and conclusion. — When the magnetic flux 
in the SQUID loop of each charge qubit is set to <I>e = 
$0/2, then hx = and Jy = 0. Also, we assume 
that Ec > Jx,Jz- When $e = <&o/2, one can shift 
the gate voltage, at the ith lattice point, far away from 
the usual working point rig ~ ^ of the Kitaev lattice 
for a period of time r. This yields a local z-type ro- 
tation on the ith qubit: i?f(0) — eyip[—i{hzT /h)a^] = 
exp(-j6'crf/2). When 6 = 2hzT/h = tt (where rig > \), 
i?f (tt) = — icrf , so the erf operation on the ith qubit is 
given by = e*'^/^i?|(7r), while half of the rotation is 
(erf) 5 — e*'^/'*i?f (7r/2). The corresponding inverse ro- 
tations can be achieved by shifting the gate voltage to 
Ua < \. Similarly, when Ug = ^ and $e = 0, one has 

0, hx = -Ej, and Jx = 
that Ej » \Jy\,Jz. When Ug — ^, 
flux in the SQUID loop of the ith qubit for a time r (the 
working point of this Kitaev lattice is usually at < < 
$0/2). This produces a local x-type rotation on the ith. 
qubit: i?f (6*) = e:>q;>[i(5EjT /n)ai] = exp(i6'crf /2), where 
5Ej — Ej — i£^j(<i>e). The af rotation on the ith qubit 
is erf = e-'''/^Rf{n), where 2SEjT/h = n. With both erf 
and erf rotations available for the ith qubit, the af ro- 
tation is given by — e^^^^-^afaf . Therefore, one can 
construct the operations (jf and af for generating topo- 
logical excitations by using the single-qubit rotations erf 
and erf. 

In order to obtain accurate z- and x-type single-qubit 
rotations, we assume that Ec and Ej are much larger 
than the inter-qubit coupling. Actually, this strict con- 
dition can be loosened for realistic systems. As shown 
in [Si J accurate effective single-qubit rotations can still 
be achieved using techniques from nuclear magnetic res- 
onance when the inter-qubit coupling is small (instead of 



Lg kJilililCLllJ', WllCll ILg ^ 

hx = 0, hx = —E,j, and Jx = 0. Here we now assume 

- ^ we also switch off the 



much smaller than Ec and Ej). In the vortex- and bond- 
state cases studied here, to reveal the braiding statistics 
of topological excitations, one should verify the occur- 
rence of the excited state [^2) at the crossing point of 
paths Ph and Py. This needs to distinguish [^2) from 
the ground state. Here the difference between IZ2) and 
the ground state is the phase flip, induced by the single- 
qubit rotation erf , on the charge state of the z-bond at 
the crossing point. This phase flip of the charge state 
could be measured, e.g., using state tomograph y |14j . 

In the region J'^ > J^,|Jy|, it is shown [15l| that 
when acting on the ground state, the spin-pair oper- 
ations af and generally create both vortex states 
(Abelian anyons) and fermionic excitations. However, 
the weight of the anyons is dominant in the low-energy 
subspace if the number m of the spin-pair operations 
for braiding anyons is not large and ^ , | Jy | . Let 
■^x — l-^yl ~ For example, when m = 25, 7 = J/J'z 
can be 7 - 0.2 [15]. 

In conclusion, we have proposed an approach to real- 
ize the Kitaev model on a honeycomb lattice using su- 
perconducting quantum circuits. Two particular cases 
are studied to demonstrate the topological states and 
the braiding statistics. Our approach provides an experi- 
mentally realizable many-body system for demonstrating 
exotic properties of these topological phases. 
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